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1 ABSTRACT 

d 

If the number of runs in a (mixed-level) orthogonal array of strength 2 is specified, what numbers 

of levels and factors are possible? The collection of possible sets of parameters for orthogonal arrays 

Qs, ' with N runs has a natural lattice structure, induced by the "expansive replacement" construction 

ON ' 

{SJ ■ method. In particular the dual atoms in this lattice are the most important parameter sets, since 

in , 

any other parameter set for an TV-run orthogonal array can be constructed from them. To get 
a sense for the number of dual atoms, and to begin to understand the lattice as a function of 



N, we investigate the height and the size of the lattice. It is shown that the height is at most 



43 

[c(N — 1)J , where c = 1.4039 . . ., and that there is an infinite sequence of values of N for which 
j> , this bound is attained. On the other hand, the number of nodes in the lattice is bounded above by 

•i-H . 

a superpolynomial function of N (and superpolynomial growth does occur for certain sequences of 
values of N). Using a new construction based on "mixed spreads", all parameter sets with 64 runs 
are determined. Four of these 64-run orthogonal arrays appear to be new. 
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1. Introduction 



Although mixed- level (or asymmetrical) orthogonal arrays have been the subject of a number of 
papers in recent years (see Chapter 9 of Hedayat, Sloane and Stufken, 1999, for references), it seems 
fair to say that we know much less about them than about fixed-level orthogonal arrays (in which 
all factors have the same number of levels). For example, there is no analogue for mixed orthogonal 
arrays of one of the most powerful construction methods for fixed-level arrays, that based on linear 
codes (see Chapters 4 and 5 of Hedayat, Sloane and Stufken, 1999). 

Again, there are many instances where the linear programming bound for fixed-level orthogonal 
arrays gives the correct answer for the minimal number of runs needed for a specified number of 
factors. There is a linear programming bound for mixed arrays (Sloane and Stufken, 1996), but it 
is less effective than in the fixed-level case — it ignores too much of the combinatorial nature of the 
problem (especially when the levels involve more than one prime number), and, though generally 
stronger than the Rao bound, does not give correct answers as often as in the fixed-level case. 

A mixed orthogonal array OA(N, s^s^ 2 • • • Sy v , t) is an array of size N x k, where k = k\ + 
&2 + • • • + k v is the total number of factors, in which the first k\ columns have symbols from 
{0, 1, . . . , s\ — 1}, the next k?, columns have symbols from {0, 1, . . . , S2 — 1}, and so on, with the 
property that in any N x t subarray every possible i-tuple of symbols occurs an equal number of 
times as a row. We usually assume 2 < s\ < S2 < ■ ■ ■ and all k{ > 1. Except in Section 5, only 
arrays of strength 2 will be considered, and we will usually omit t from the symbol for the array. 

We refer to (N, s^ 1 s^ 2 . . .) as the parameter set for an OA(N, s^s^ 2 . . .). We also allow the 
parameter set (N, l 1 ), corresponding to the trivial array consisting of a single column of N 0's. 
In this paper we consider the question: if iV is specified, how many different parameter sets are 
possible? 

Given an array A = OA^s^s^ ...), other N -run arrays can be obtained from it by the expan- 
sive replacement method. Let S be one of the Sj occurring in A, and suppose B is an OA{S, 
The expansive replacement method replaces a single column of A at S levels by the rows of B. 
For example, if A = OA(16, 2 3 4 4 ) and B = OA(4, 2 3 ), we obtain an OA(W, 2 6 4 3 ). If B is a trivial 
array OA(S, l 1 ), we are simply deleting one of the S'-level factors from A. E.g. taking S = 2, 
an OA(2A, 2 20 4 1 ) trivially produces an OA(24, 2 19 4 1 ). The expansive replacement method also in- 
cludes replacing a factor at s levels by a factor at s' levels, if s' divides s. For further details about 
the expansive replacement method see Chapter 9 of Hedayat, Sloane and Stufken, 1999. 

Let A and B be parameter sets for orthogonal arrays with iV runs. We say that B is dominated 
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by A if an orthogonal array with parameter set B can be obtained from an orthogonal array with 
parameter set A by a sequence of expansive replacements. 

Using "dominance" as the relation, the parameter sets for orthogonal arrays with N runs form 
a partially ordered set, which we denote by A^r (Hedayat, Sloane and Stufken, 1999 p. 335). 

Atv has a unique maximal element (NjN 1 ) (corresponding to the trivial array with one factor 
at N levels) and a unique minimal element (N, l 1 ). It is straightforward to verify that meet (A) 
and join (V) are well-defined for this relation (we omit the proof), so Ajy is in fact a lattice (cf. 
Welsh, 1976; Trotter, 1995). 

If an OA(N, s^ 1 ^ 2 . . .) exists, then necessarily we must have: 

(CI) Si divides N, for all i, 
(C2) sf divides N, if > 2, 
(C3) S{Sj divides N, if i ^ j, 
(C4) the Rao bound holds: 

N-l> fci(ai-l) + k 2 {s 2 + , (1) 
(C5) the linear programming bound holds (see Sloane and Stufken, 1996). 

These conditions are certainly not sufficient for an array to exist, and it appears to be difficult 
to test if an orthogonal array does exist with a putative parameter set satisfying (C1)-(C5). A 
further difficulty is that in order to construct Ajv it is necessary to know A^ for all proper divisors 
d of N. 

To avoid these difficulties we define a second lattice, the idealized lattice A' N : this has as nodes 
all putative parameter sets satisfying (CI) to (C4), with the dominance relation as before, except 
that in the expansive replacement method we may now make use of any of the nodes of any A^ for 
d dividing N . 

Constructing A'^ is much easier than constructing Ajy, since essentially all we need to do is 
enumerate the solutions to (||). Of course A^v is a sublattice of A^. 

To avoid having to repeat the adjective "putative", from now on we will use "parameter set" 
to mean any symbol (N, s^ 1 s 2 2 ■ ■ ■) satisfying conditions (CI) to (C4). The parameter sets are 
precisely the nodes of A^y. If a parameter set is also a node of Ajv then it is implied that an 
OA(N, s^ 1 ^ 2 • • •) does exist, i.e. that the parameter set is realized by an orthogonal array. 

It is convenient to represent A^v and A^y by their Hasse diagrams (cf. Welsh, 1976, p. 45). 
These diagrams are drawn "from the bottom up", with (A r , l 1 ) as the root node at the bottom 
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(Figure [l] shows A12 and A' 12 ). The height of a parameter set is the number of edges in the longest 
path from that node to the root. A node of height i appears on the zth level of the diagram. The 
height of the maximal element (N, N 1 ) will be denoted by ht(N). 



12 1 12 1 




Figure 1: (Left) Lattice A12, showing all (true) parameter sets for 12-run orthogonal arrays. There 
are 23 nodes, four dual atoms (circled) and the height is 12. (Right) Idealized lattice A' 12 , showing 
all (putative) parameter sets satisfying conditions (C1)-(C4). No arrays exist for the nine nodes 
marked x . There are 32 nodes, four dual atoms (circled) and again the height is 12. 

The atoms in A^r (those nodes just above the root) are precisely the parameter sets (A^p 1 ) for 
the primes p dividing N. 

The dual atoms in Ajv (those nodes just below the maximal element) are especially interesting, 
since they dominate all other parameter sets. 

We can now state our main results. 

Theorem 1. (i) For all N , 

ht{N) < [c{N - 1)J , (2) 

where 

00 

c = E^3T = L4039 --- • (3) 
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(%%) IfN = 2 2m (m > 0) then ht(N) = [c(N - 1)J . 

Let T(N) (resp. T'(N)) denote the total number of nodes in Ayv (resp. A^). 
Theorem 2. If N = 2 n , 

^(log 2 iV) 2 (l + o(l)) < log 2 T(JV) < ^(log 2 iV) 2 (l + o(l)). (4) 

Theorem 3. There is a constant c\ such that for all N, 

lnlnr(AT) < Cl (l + o(l)). (5) 

Remarks, (i) The bounds in @ and @ also apply to T'(N). 

fii) Theorem | shows that when A" = 2™, T(N) grows very roughly like N alog2N , for some 
constant a between | and |. This is a "superpolynomial" function of N, meaning that it grows 
faster than any polynomial in N. 

(iii) It appears (although we have not proved this) that the upper bound in (@) can be achieved 

by taking A^ to be a certain product of powers of the first m primes, where m is about 

1 In A" 
Ye In In N 

(see Section 7). In other words, it appears that there is an infinite sequence of values of A^ for 
which T(N) grows very roughly like 

exp(A^ C2 / lnln7V ) , 

where C2 is a constant. This is again a superpolynomial function of N, and is now close to being 
an exponential function, since lnlnAf grows slowly. 

The above discussion has shown that there is an infinite sequence of values of A^ for which the 
number of nodes in A^r grows superpolynomially, while the height of A^v grows at most linearly. It 
follows that the size of the largest antichain must also grow superpolynomially. The data in Table 
H suggest the following conjecture. 

Conjecture. There is an infinite sequence of values of N for which the number of dual atoms 
grows superpolynomially in N. 

In fact it seems likely that if A^ = 2 n , a lower bound of the form in (||) (possibly with a different 
constant) applies to the logarithm of the number of dual atoms, and that for some sequence of 
values of A^ a lower bound similar to the upper bound on the right-hand side of (|B|) will hold. 
However, at present these are only conjectures. 
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In order to construct the orthogonal arrays needed to establish the lower bound in Theorem || 
we make use of what we call "mixed spreads", generalizing the notions of "spread" and "partial 
spread" from projective geometry. Arrays that can be constructed in this way we call "geometric" . 
Many familiar examples of orthogonal arrays, for example arrays constructed from linear codes, 
are geometric. The construction is not restricted to strength 2 (and is one of the few general 
constructions we know of for mixed arrays of strength greater than 2). The construction will be 
described in Section 5. 

In Section 6 we use this construction to determine the lattice Ag4, and in doing so we find tight 
arrays with parameter sets 

(64, 2 5 4 17 8 1 ), (64, 4 14 8 3 ), (64, 2 5 4 10 8 4 ), (64, 4 7 8 6 ), 

which appear to be new. 

When studying parameter sets of putative orthogonal arrays with iV runs, it is convenient to 
be able to say that if the number of degrees of freedom of the parameter set (N, s^s^ 2 • • •), that is, 

fci(ai-l) + *2(s2-l) + "- , (6) 

is small compared with N — 1, then an orthogonal array certainly exists. 

To make this precise, we define the threshold function B(N) to be the maximum number b such 
that every parameter set (satisfying conditions (CI) to (C4)) with at most b degrees of freedom is 
realized by an orthogonal array, but some parameter set (again satisfying (CI) to (C4)) with 6+1 
degrees of freedom is not realized. If every parameter set satisfying (CI) to (C4) is realized, we set 
B(N) =N-1. 

Figure [l] shows that -6(12) = 6, since there is no OA(12, 2 5 3 1 ), but every parameter set with at 
most 6 degrees of freedom is realized. 

We are not aware of any earlier investigations of B(N). 

Theorem 4. If N is a power of a prime then 

iV 3/4 < B(N) . 

In words, if the number of degrees of freedom in the parameter set does not exceed iV 3 / 4 , then 
an orthogonal array exists. This is certainly weak, but is enough to establish the lower bound of 
Theorem ||. It would be nice to have more precise estimates for B(N). 

A final remark. We could have considered the partially ordered set whose nodes are all the 
inequivalent orthogonal arrays with iV runs, rather than just their parameter sets. However, the 

6 



number of nodes then becomes unmanageably large, even for small values of iV (furthermore, it 
appears that "meet" and "join" are no longer well-defined, and so in general this partially ordered 
set would not be a lattice). 

Consider N = 28, for example. Using Kimura's (1994a, 1994b) enumeration of the Hadamard 
matrices of order 28, we have calculated^ that there are precisely 7570 inequivalent 0^4(28, 2 27 )'s. 
This would be merely a lower bound on the number of dual atoms. On the other hand we know 
(see Table |l]) that A28 has precisely four dual atoms, between 47 and 55 nodes, and height 28. 

(a) 



l 1 



A, 




Figure 2: (a) A p and (b) Aq. 



2. Examples of the lattices A at and A', 



N 



There are a few general cases when we can describe A^r explicitly (and for which A^ is the same 
as Ajy). 

If N = p is a prime then A^v = A^ has two nodes, one dual atom and height 1, as shown in 
Fig. ||(a) (dual atoms are circled). 

If N = pq is the product of two distinct primes, Ajy = A'^ has five nodes, one dual atom and 
height 3. Aq is shown in Fig. ||(b). 

More generally, if N is the product of u > 2 distinct primes, it is not difficult to show that 
Ajy = A' N has 2" _1 — 1 dual atoms, height 2u — 1, and u +i nodes, where 



{ft, ft, /%,...} = {1,2,5,15,52,...} 



are the Bell numbers (see Sequence A11C of Sloane, 1999). Figure y shows A30, illustrating the case 
u = 3. 

Ii N = p 2 is the square of a prime, A^r = A^ has a single dual atom, OA(p 2 ,p p+1 ), has height 

p + 2 and contains p + 3 nodes. If N = p 3 (p prime), A^r = A^ also has a single dual atom, 
1 Using the method described on page 165 of Hedayat, Sloane and Stufken (1999). 



30 1 




1 1 



Figure 3: The lattice A30: there are (3$ = 15 nodes, three dual atoms (circled) and the height is 5. 

OA(p 3 , (p) p2 (p 2 ) 1 ), has height p 2 + p + 3 and contains 2p 2 + p + 4 nodes. As is shown in Fig. [|. 
If N = j> 4 (p prime), An = A' N has two dual atoms, OA(p 4 , (p 2 ) p2+1 ) and OA(p^, (p) pS (p 3 ) 1 ), has 
height p 3 + 2p 2 + p + 3 and contains 

-(p 5 + p 4 + 5p 3 + hp 2 + 2p+ 10) 

nodes. 

For all values of N mentioned so far in this section, the threshold function B(N) = N — 1. 

If N is not of one of the above forms then it seems necessary to consider each case individually. 
Table [l] summarizes the properties of Ayv for some small values of N. Here A(N) denotes the 
number of dual atoms in A^v- Most of the entries in this table can be deduced from the tables in 
Chapter 12 of Hedayat, Sloane and Stufken (1999), except for N = 32 and 64, which are discussed 
in Section 6. 

N = 24 is the smallest case when we do not know the complete lattice Ayv- In A24 the maximal 
value of k that occurs in each of the parameter sets 2 fc 3 1 , 2 fc 3 1 4 1 , 2 fc 4 1 6 1 and 2 fc 6 1 is presently 
unknown. For 2 fc 3 1 , for example, it is known only that an OA(2A, 2 16 3 1 ) exists, and an OA(2A, 2 21 3 1 ) 
is impossible by the linear programming bound. The number of dual atoms is in the range 4 to 7. 
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2 7 



2- 

2 1 
l 1 



Figure 4: The lattice As: 14 nodes, one dual atom, height 9. 



It is possible to show that the height of A24 is 25, however: no chain can be longer than 



24 1 _ 2 20 4 1 _ 2 23 _ 2 22 _ 2 21 



2 1 - l 1 . 



We also do not know for N = 28, 36, . . .. 

The four sequences in Table [l] are Sequences |A39927| , |A3993q , |A39931| and |A48893| of Sloane 
(1999). The entries in that database will be updated as further values are determined. 

We end this section with a remark about the nature of Ajy as a lattice. Since not all maximal 
chains between two nodes need have the same length (see Fig. [l]), A^r does not in general satisfy the 
Jordan-Dedekind chain condition (cf. Welsh, 1976). It follows that A^r is in general not distributive, 
not semimodular, nor is it the lattice of a matroid. 

3. The maximum height of A^r 

In this section we give the proof of Theorem ||. 

Let a denote a specification s^s^ 2 ... of factors at various levels, leaving the number of runs 
unspecified. Given a, there is a smallest number of runs, Nq say, for which an OA(Nq,(j) exists. 
Let h be the height of the parameter set (No, a) in Atv - Then if the parameter set (N,a) occurs 
in any other lattice A^r, it also has height h. (E.g. the specification a = 6 1 has height 3 in each of 
Figs. 1, 2(b) and 3.) We may therefore define ht(a) to be h, independently of the number of runs. 



9 



Table 1: For the lattice A^v of parameter sets of orthogonal arrays with N runs, the table gives the 
number of dual atoms A(N), the height ht(N), the total number of nodes T(N) and the threshold 
function B(N). 



N 


A(N) 


ht(N) 


T(N) 


B(P 


1 








1 





2 


1 


1 


2 


1 


3 


1 


1 


2 


2 


4 


1 


4 


5 


3 


5 


1 


1 


2 


4 


6 


1 


3 


5 


5 


7 


1 


1 


2 


6 


8 


1 


9 


14 


7 


9 


1 


5 


6 


8 


10 


1 


3 


5 


9 


11 


1 


1 


2 


10 


12 


4 


12 


23 


6 


13 


1 


1 


2 


12 


14 


1 


3 


5 


13 


15 


1 


3 


5 


14 


16 


2 


21 


61 


15 


17 


1 


1 


2 


16 


IS 


2 


10 


26 


15 



The height of A^v, ht(N), as defined in Section 1 coincides with ht(a) when a = N . The height 
function has the following additivity property. 

Lemma 5. 

M^s* 2 ...) = k 1 ht(s 1 ) + k 2 ht(s 2 ) + --- . (7) 

Proof. If there is a single factor on the left-hand side, say s± = s, k\ = 1, then just says that 
ht{s l ) = ht(s), repeating the assertion made above. Otherwise, more than one factor occurs in the 
specification a on the left-hand side of (0). Suppose the parameter set (N,a) occurs as a node in 
some Ajy. The portion of A^ consisting of this node and all nodes dominated by it is the product 
in an obvious sense of k\ copies of A S1 , k 2 copies of A S2 , etc. The height of (N, a) is the sum of the 
heights of all these sublattices, and (0) follows. ■ 

Lemma |5] reduces the calculation of ht^^s 1 ^ 2 • • •) to the calculation of the values of ht(si), 
ht(s 2 ), .... To determine the latter we must consider exactly which orthogonal arrays exist with a 
specified number of runs. For this we can apply the following lemma. 
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Table 1 (cont.) 



N 


A(N) 


ht(N) 


T(N) 


5(JV) 


19 


1 


1 


2 


18 


20 


4 


20 


35 


11 


21 


1 


3 


5 


20 


22 


1 


3 


5 


21 


23 


1 


1 


2 


22 


24 


4-7 


25 


119 - 133 


18-22 


25 


1 


7 


8 


24 


26 


1 


3 


5 


25 


27 


1 


15 


25 


26 


98 


A 
-t 






1 o 


29 


1 


1 


2 


28 


30 


3 


5 


15 


29 


31 


1 


1 


2 


30 


32 


2 


42 


320 


29 


33 


1 


3 


5 


32 


34 


1 


3 


5 


33 


35 


1 


3 


5 


34 


64 


7 


86 


3037 


57 



Lemma 6. 

ht(N) = 1 + max^2kiht(si) , (8) 

i 

where the maximum is taken over all parameter sets (N, s^ 1 ^ 2 • • •) 7^ (-^> N 1 ) for which an orthog- 
onal array exists. 

Proof. The height of A^r is one more than the maximal height among the dual atoms. @ follows 
by applying Lemma [| to the parameter set of such a dual atom. ■ 

We can now use linear programming to obtain an upper bound on ht(N), by maximizing 

l + Y, k iht( Si ) (9) 

i 

over all choices of s\, k±, S2, &2, • • • that satisfy (CI) to (C4). 
We first consider the case when N = 2 n for some n. 

The case iV = 64 will illustrate the method. If there is a factor 32 1 then linear programming 
shows that (|9|) is maximized by 2 32 32 1 , giving height 75. If there is a factor 16 1 then there is 
a unique parameter set that maximizes (|9|), 4 16 16 1 , giving height 86. Otherwise, if only factors 
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2 fcl 4 fc2 8 fca occur, the height does not exceed 85. Since an OA(64, 4 16 16 1 ) exists, we conclude that 
ht{QA) = 86. We will return to the case N = 64 in Section 6. 

In this way we obtain the values of ht(N), N = 2 n shown in Table |^. 

Table 2: Height of for N = 2 n . 



71 


1 


2 


3 


4 


5 


N = 2 n 


2 


4 


8 


16 


32 


ht(N) 


1 


4 


9 


21 


42 


ht(N)/(N - 1) 


1 


1.3333... 


. 1.2857... 


1.4 


1.3548. 


n 


(3 


7 


8 


9 


10 


N = 2 n 


64 


128 


256 


512 


1024 


ht{N) 


86 


171 


358 


715 


1431 


ht(N)/(N - 1) 


1.3650. . 


. 1.3465.. 


, 1.4039... 


1.3992... 


. 1.3988. 



Consider the general problem of maximizing @ for N = 2 n . Comparing (|T]) and (|9|), we see 
that an s-level factor contributes ht(s) to the height but uses up s — 1 degrees of freedom. If we 
ignore the constraints of integrality then the value of the expression in @ would be maximized by 
a term s k where s is chosen to maximize ht(s)/(s — 1). This suggests that we should investigate 
this quantity in order to prove Theorem |l[ The data in Table ^ suggest that the ratio ht(s) /(s — 1) 
is maximized if s is of the form 2 2 ' and is as large as possible. That this is indeed so is established 
by the next three lemmas. 

Lemma 7. Given positive real numbers a r , (3 r (r = 1, . . . ,m), 7, the maximal value of 

in 

^^^^^ OLf X y 



r=l 



subject to the constraints 



PrX r = 7 , 



r=l 



x r > (r = 1, . . . , m) 



IS 



7 max — 

r=l...m p r 



We omit the straightforward proof. 
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Lemma 8. Given positive real numbers a r , j3 r (r = 1, . . . , m + n), 7 with 7 < 7, the maximal 
value of 

m+n 

a r x r (10) 



r=l 



subject to the constraints 



m+n 

^2 (3 r X r = 7 
r=l 



m+n 



r=m+l 

x r > (r = 1, . . . , m + n) 
max < 7 max — , (7 — 7 J max — — h 7 

r=l...m p r r=l...m p r r= 

The maximum is given by the first expression if and only if 



max — 

-m+l...m+n p r 



(11) 

(12) 
(13) 

(14) 



max — > max — 
r=l...m p r r=m+l...m+n p r 



Proof. Let a, < a < 7', denote the value of the left-hand side of (|l^). Then by Lemma ^ the 
maximal value of the sum in ( |To|) is equal to 



(7 — a) max — + a max — . 
r=l...m p r r=m+\...m+n jj r 

This is a linear function of a and so its maximal value is taken at one of the two endpoints, leading 
to©- ■ 

We can now give an upper bound on the height ht(N) for N = 2 n . Let 

ht{N) 



Pi. 



2™ - 1 



N = T 



Lemma 9. 



where c is the constant 



p n < c for all n 



00 ^ 

o2^ _ 1 
i=0 z 1 



1.4039. 
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Proof. We use the induction hypothesis that, for m > 0, 

p 2 m > p x whenever 1 < x < 2 m+1 , x + 2 m . (15) 

This is trivially true when m = 0. We first compute p 2 m. Prom the Rao-Hamming construction, 
for TV = 2 2m , an OA(N, (y/N) VJ!+1 ) always exists, and we find that 

1 

On the other hand, we obtain an upper bound on p 2 m from the linear program: choose nonnegative 
integers k\ , k 2 , . . . so as to maximize 

2 m -l 

1 + ht(2 r )k r 

r=l 

subject to the constraint 

2 m -l 

(2 r - l)K = 2 2m - 1 . 

From Lemmas ^ and |^ we have 



r=l 



P2m < ^ r - T+r= max_ /r 



+ p 2 m-l . 



2 2m - 1 
Since the two bounds agree, 

p 2m = p 2m -i + 22m 1 _ - . (16) 
We now complete the proof of the induction step. For 1 < x < 2 m , we have 

p 2 m > p 2 m-l > p x , 

as required. Suppose 2 m < x < 2 m+1 . Then p x is upper-bounded by the solution to the linear 
program: maximize 

(l + J2ht(2 r )k r ) /(2 X -1) 



r=l 



subject to the constraints 



x-l 



^2(2 r - l)k r = 2 X - 1 , 

r=l 
x-l 

(2 r - l)k r < 2 x ~ l - 1 . 



r=2" 
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(The second constraint is implied by the requirement that there can be at most one factor which 
has more levels than the square root of the number of runs.) By Lemma || and induction on x, we 
obtain 

1 



VS— 1 , fr,X— 1 

i fJ^TYl—l 

1 + 2*-i - 2 2 ' 



< ^— T (l + 2-V 2m -i + (2 a - 1 -l)p 2m ) 



(2 2m - 1)(2* - 1) 

< P2 m ■ 

This completes the induction step. 

To complete the proof of Lemma pi by the induction hypothesis it suffices to prove that p2 m < c 



for all m. But from (16) it follows that 

^ 2 2 ' - 1 



p 2m = ^ — < c , 

i=0 

and that /?2 m — > c as m — > oo. ■ 

If N is not a power of 2 then similar arguments show that the height is (considerably) less than 
cN. This completes the proof of Theorem ||. 

4. Upper bounds on the number of parameter sets 

In this section we establish the upper bounds in Theorems |2] and |3|. We will bound T'(N), the 
number of nodes in A' N . Since Ajy is a sublattice of A^y, this is also an upper bound on the number 
of nodes in An. Suppose first that N = 2 2r . 

We start by considering parameter sets (N, 2 kl 4 k2 8 ki . . . (2 r ) kr ), containing no level exceeding 
VN. From (0), 

N - 1 > h + 3k 2 + 7k 3 + • • • + (2 r - l)kr . (17) 

Let 7 be the number of nonnegative integer solutions (fci, fe, . . . , fcjv) to this inequality. Then 
"f/(N — l) r is the Riemann sum approximating the volume of the simplex bounded by the hyper- 
planes 

1 > Xi + 3x 2 + 7x 3 + • • • + (2 r - l)x r , 
X\ > 0, X2 > 0, . . . , x r > 

in W . For large this yields 

1 = ^^- (l + o(l)). (18) 

r! n(2*-l) 

i=i 
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The product in the denominator approaches c^2 r ^ r+l ^ 2 , as r — > oo, where C3 = 0.2887 .... 

Now suppose the parameter set contains a factor at 2 % levels, where r + l<i<2r — 1. There 
can be at most one such factor, and the number of such parameter sets in each case is at most 7. 
The total number of parameter sets is therefore at most r^, and setting r = ^ log 2 N we find that 

log 2 (r 7 )<^(log 2 A0 2 (l + o(l)). 

This establishes the upper bound in Theorem || for N = 2 2r . It also implies the upper bound for 
N = 2 2r+1 , after noting that T'(N) < T'(2N). 

We now give a sketch of the proof of Theorem ||, omitting many tedious details. To simplify 
the analysis we will neglect terms on the right-hand side of (||) that correspond to factors with a 
level greater than y~N. Suppose first that N is a large number of the form 2 2ai 3 2a2 . Pretending 
for the moment that a\ and a 2 are allowed to be real numbers, not just integers, we may consider 
what choice of a\ and a 2 maximizes the number of solutions to ([!]) for a given value of N. The 
number of terms on the right-hand side of (|l|) is now (a\ + l)(a 2 + 1) — 1. The arguments used to 
establish the upper bound of Theorem || show that the number of solutions to (|l]) is maximized if 
2 2<Jl is approximately equal to 3 2ct2 . 

Now suppose that N is of the form 

Vl ai V2 2 ---P 2 « m , (19) 

where p\ = 2, p 2 = 3, . . . are the first m primes. We find that the number of solutions to (0) is 
maximized when the numbe 
p 2 ^ = jv£(i+°(i)), i.e. that 



maximized when the numbers p 2ai are all approximately equal, and we will therefore assume that 



1 In 

ai = — (l + o(l)), i = l,. ..,m. 

2m In pi 

The Rao bound contains a term for every possible level 

b — P\ P2 Pm 

in which < % v < a u , 1 < v < m, where not all the i u are equal to 0. The number of such terms is 

1 0nA0 m 

5 := ( Ql + l)(o2 + 1) ■ • • (a m + 1) - 1 = V m ' (1 + o(l)) . (20) 

n 

3=1 

The product of the coefficients of all the terms on the right-hand side of the Rao bound is 

C . = p ala 2 ...a m /2 p a ia la 3 ...a m /2 _ _ _ ^ + _ 
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This implies 

InC = -InJV(l + o(l)) . 

Again using 7 to denote the number of solutions to the Rao inequality, we have 

N 5 

7 = ^(1 + 0(1)), 

hence 

3 (lniV) m+1 1 (lniV) m 

hi 7 = —- ; — ; : — r- — tttt — 7 ; — ~, ; — r- — tttt(^ In In N — m In m) 

4 (2m) m (lnm) m ( 1 +°( 1 )) (2m) m (lnm) m d+°(i)) v ' 

+ smaller terms . 
This expression is maximized if we take 

1 In AT , 

m= 2^hThTiv (1 + o(1)) ' 

and then we find that the leading term in the expression for In 7 is 

3 1 1 

— In iV N 2e In In JV 

4 



1 IniV 

In In 7 < 



2elnlniV 



We conclude that 

which establishes Theorem [3]. 

5. Geometric orthogonal arrays 

We consider subspaces V of the vector space GF(q) n over GF(q), where q is a power of a prime. 
By the dimension of V, dimV, we mean the vector space dimension over GF(q) (rather than the 
projective dimension, which is one less). The following notion was suggested by the notions of 
spread and partial spread in projective geometry (cf. Thas, 1995). 

Definition. A mixed spread of strength t is a collection V = {V±, V2, ■ ■ ■ , V^] of subspaces of 
GF{q) n such that for all choices of r < t indices i\, Z2, . • • , i T (with 1 < i\ < ii < ■ ■ ■ < i T < k) the 
dimension of the span of , . . . , V% T is equal to dim + • • • + dim Vi T . 

An equivalent condition is that the span of , . . . , Vi T is the direct sum © • • • © Vi T for all 
choices of r < t indices i±, . . . , i T with 1 < i\ < ■ ■ ■ < i T < k. 

Any collection V of subspaces has strength 1. V has strength 2 if and only if every pair Vi, 
Vj E V, i 7^ j, intersect just in the zero vector. V has strength 3 if and only if it has strength 2 and 
for any triple of distinct subspaces each one meets the span of the other two just in the zero vector. 

17 



If V is a c?-dimensional subspace of GF(q) n we denote by V* the dual space, the space of linear 
functionals on V (see for example Hoffman and Kunze, 1961), and we fix a labeling /o, /i, . . . , fqd-i 
for the elements of V* . 

Given a mixed spread of strength t, V = {Vi, V2, . . . , Vk}, where the Vi are subspaces of GF(q) n , 
we obtain an orthogonal array OA(V) with q n runs as follows. The columns of the array are labeled 
Vi, V2, • • • , Vfc and the rows are labeled by the linear functionals / £ (GF(q) n )* . If / restricted to 
Vi, f\vi, is the jth linear functional in V*, the (/, VJ)-th entry in the array is j. The symbols in 
column i are therefore taken from {0, 1, ... , q dimVi — 1}. 

We will say that an orthogonal array constructed in this way is geometric. 

Theorem 10. The orthogonal array OA(V) has strength t if and only if the mixed spread V has 
strength t. 

Proof. Suppose OA(V) has strength t. Consider for example the first t columns. In the projection 
of the array onto these columns we see 

t ' 
J q dimVi =9 £i 

i=l 

different t-tuples of symbols. Since these depend only on the restrictions of the / G {GF{q) n )* to 
the span of V±, . . . , Vt, the dimension of that space must be at least the sum of the dimensions of 
V±, . . . ,Vt, and clearly it cannot have a higher dimension. So V is a mixed spread of strength t. 
Conversely, suppose V is a mixed spread of strength t. We can write 

GF(q) n = V l @---@V t @X 

where X is the complementary space to the Vi- Since the dual of a direct sum is canonically 
isomorphic to the direct sum of the duals, we immediately find that as we run through the linear 
functionals on GF(q) n , every tuple . . . /|y t , f\x) of restrictions occurs precisely once. Ignoring 

the last component, we see that every tuple (/|vi, . . . f\v t ) occurs precisely \X\ times. Hence OA(V) 
has strength t. m 

Lemma 11. Any geometric array of strength 2 can always be extended to a tight array (i.e. one 
meeting the Rao bound) by adding q-level factors. 

Proof. We simply group any unused points into 1-dimensional subspaces. ■ 
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Examples, (i) The 1-dimensional subspaces of GF(q) n form a mixed spread of strength 2. The 
corresponding array is the familiar 

OA(q n ,q k ), k = (<f - l)/(q - 1) , 

of the Rao-Hamming construction. 

(ii) More generally, a classical a-spread in PG(b, q) is a mixed spread of strength 2 in our sense. 
This is a set of subspaces of PG(b, q) of projective dimension a which partitions PG(b, q) (Thas, 



1995), and exists if and only if a + 1 divides 6+1. From Theorem 1C we obtain an 



OA{q b+ \ (q a+1 ) k ), k = (q b+1 - l)/(q a+1 - 1) , 

which of course is also given by the Rao-Hamming construction. 

We could also have obtained example (ii) directly from example (i), by remarking that a mixed 
spread of strength t over GF(q), q = p@, is also a mixed spread of strength t over GF(q > ), q' = p a , 
if q' divides q. The dimensions of the subspaces are multiplied by (3 /a. 

(iii) Provided a > 6/2, there exists a mixed spread of strength 2 in GF(q) b consisting of a single 
subspace GF(q) a and a partitioning of the remaining points into q a subspaces GF(q) b ~ a . This 
can be proved directly, or alternatively is equivalent to Lemma 2.1 of Eisfeld, Storme and Sziklai 



(1999). From Theorem 10 we obtain a geometric 

OA(q b ,(q b - a y a (q a ) 1 ) 

whenever a > 6/2. Orthogonal arrays with these parameters were already known from the dif- 
ference scheme construction (Hedayat, Sloane and Stufken, 1999, Example 9.19), but the present 
construction also shows that they are geometric. 

(iv) The classical "partial a-spread" constructed in Lemma 2.2 of Eisfeld, Storme and Sziklai 
(1999) translates in our language into a mixed spread of strength 2 consisting of k 6-dimensional 
subspaces (6 > 2) of GF(q) n , where n = ib + r, < r < 6, and 

a ib - 1 
k = q r ? 1 —-q r + l. 
q° — 1 

This produces a geometric OA(q n , (q b ) k ) (again arrays with these parameters were known from the 



difference scheme construction), which by Lemma 11 can be extended to a tight 



OA{q n ,{q X Y{<fr) , (21) 

where I = q b (q r — l)/(q — 1). The orthogonal arrays constructed by Wu (1989) are a special case 
of (pH), but in general these arrays may be new. 
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(v) Generalizing examples (i) and (ii), any orthogonal array formed from the codewords of a 
projective linear code (one for which the columns of a generator matrix are nonzero and projectively 
distinct) is geometric. 

(vi) The 0^4(256, 2 16 ) of strength 5 formed from the Nordstrom-Robinson code (see Hedayat, 
Sloane and Stufken, 1999, Section 5.10) is not geometric, and no geometric 0^4(256, 2 16 ) of strength 
5 exists. 

We shall see other examples in Section 6. 

Remarks. An unmixed geometric orthogonal array is always linear, in the sense of Hedayat, 
Sloane and Stufken (1999), Chapter 3. In general a mixed geometric orthogonal array is additive 
but not necessarily linear^] over each of the fields involved. 

If the strength is 2, the number of degrees of freedom in the parameter set for OA(V) is equal 
to the total number of nonzero points in all the subspaces Vi. 

Finally, the following is a recipe for constructing the orthogonal array from a mixed spread 
V = {Vi, Vz, ■ ■ ■ , Vk} of subspaces of GF(q) n in the case when q is a prime. Let v± , . . . , be a 
basis for Vi, where dj = dim Vi, 1 < i < k. Let Wq, . . . , u? g n_i be the vectors of GF(q) n . Then the 
ith entry of the jth row of the orthogonal array, for 1 < i < k, < j < q n — 1, is the number 

X>^« tf- 1 . 

r=l 

(This is a number in the range {0, . . . , q di — 1}.) 
6. If the number of runs is a power of 2 

In this section we consider the case N = 2 n , n = 1,2,.... We have already discussed ht(N) in 

Section 3 (see Table ||). With the assistance of Michele Colgan, we used a computer to determine 

the number of dual atoms ^4' (iV) and the total number of nodes T'(N) in the idealized lattice 

for n < 9. The results are shown in the second and third columns of Table ||. Note in particular 

the extremely rapid growth from N = 256 to N = 512. We regard this as convincing evidence that 

when N = 2 n , A'(N) (and therefore presumably ^4(iV)) grows faster than any polynomial in N. 

As to the lattice Atv itself, for n < 4 this is covered by the results in Section 2. For N = 32 

there are precisely two parameter sets in A' 32 which do not exist, (32, 4 10 ) and (32, 2 1 4 10 ). These 

can be ruled out either by the linear programming bound or by the Bose-Bush bound (Hedayat, 
2 For the distinction between additive and linear sets in the context of coding theory see Calderbank et al. (1998). 
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Table 3: Dual atoms and total number of nodes in idealized lattice A^y (A'(N) and T'(N)) and in 
lattice Aat (A(N) and T(N)). 



N 


A'(N) 


T'(N) 


A(N) 


T(N) 


1 





1 





1 


2 


1 


2 


1 


2 


4 


1 


5 


1 


5 


8 


1 


14 


1 


14 


16 


2 


61 


2 


61 


32 


3 


322 


2 


320 


64 


11 


3058 


7 


3037 


128 


21 


33364 






256 


72 


789085 






512 


144521 


18614215 







Sloane and Stufken, 1999, Theorem 2.8). All other parameter sets in Ag 2 are realized. It follows 
that A 32 contains exactly two dual atoms, OA(32, 2 16 16 1 ) and OA(32,4 8 8 1 ). 

Before considering Ag4 we give a lemma that will be used to construct new arrays. 

Lemma 12. Suppose V\, V<i, V3 are three r- dimensional subspaces of GF(2) 2r such that ViCiVj = 
{0}; * 3- Then their union can be replaced by 2 r — 1 two-dimensional subspaces, any pair of which 
meet just in the zero vector. 

Proof. Since V\ D Vi = {0}, V\ and V2 span the space GF(2) 2r . Let m, tt2 be the associated 
projection maps from GF{2) 2r to Vi, V2 respectively. Then i± = iri\y 3 '- V3 — > V\ and 12 = 7r 2|v r 3 : 
V3 — > V2 are both isomorphisms. It follows that V3 is the set 

{v + i(v) : v in Vi}, i = i2i\ l ■ 

But then we need simply take the planes {0, v, i(v),v + i(v)} for v G V± to establish the lemma. ■ 
The lemma implies that if a geometric OA(2 2r , . . .) exists then so does the array obtained by 
replacing (2 r ) 3 in the parameter set by 4 fc , k = 2 r — 1. In particular, in a geometric 0^4(64, . . .) we 
can replace 8 3 by 4 7 . 

Theorem 13. The lattice Ag4 contains precisely seven dual atoms, with parameter sets 

(64, 2 5 4 17 8 1 ), (64, 4 14 8 3 ), (64, 2 5 4 10 8 4 ), 

(64, 4 7 8 6 ), (64, 8 9 ), (64, 4 16 16 1 ), (64, 2 32 32 x ) . (22) 
A geometric orthogonal array exists for each of these parameter sets. 
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Proof. As an intermediate step, we use mixed spreads of strength 2 to construct orthogonal 
arrays with the following parameter sets: 

(64, 4 21 ), (64, 2 5 4 17 8 1 ), (64, 2 4 4 15 8 2 ), (64, 4 14 8 3 ), 
(64, 2 5 4 10 8 4 ), (64, 2 4 4 8 8 5 ), (64, 4 7 8 6 ), (64, 2 8 4 2 8 7 ), 

(64, 2 4 4 1 8 8 ), (64, 8 9 ), (64, 4 16 16 1 ), (64, 2 32 32 1 ) . (23) 

On the other hand, linear programming shows that orthogonal arrays do not exist with parameter 
sets 

(64, 4 18 8 1 ), (64, 4 16 8 2 ), (64, 4 n 8 4 ), (64, 4 9 8 5 ), (64, 4 3 8 7 ), (64, 4 2 8 8 ) . (24) 

We then check that every parameter set with 64 runs either dominates one of (|2~| ) (and so is not 
realized), or is dominated by one of (p3|) (and is realized). Furthermore, the parameter sets in ( f22| ) 
dominate all of (f23|). 

It remains to construct the arrays mentioned in (|2^). The last two follow from Example (iii) of 
Section 5. Also (64, 8 9 ) comes from Example (i), and (64, 4 7 8 6 ) and (64, 4 14 8 3 ) follow from Lemma 
12. 

To construct an OA(64, 2 5 4 10 8 4 ) we proceed as follows. We begin by constructing an explicit 
example of an 0,4(64, 8 9 ) from Theorem [TOl, by using an extended Reed-Solomon code of length 
9, dimension 2 and minimal distance 8 over GF(8). This gives a decomposition of GF(2) G into 9 
copies of GF(2) 3 meeting only in the zero vector. These 9 subspaces are spanned by the following 
nine triples of columns: 



I I I I I I I I 

1 / A A 2 A 3 A 4 A 5 A 6 



(25) 



where 





" 





" 




" 1 





" 




' 


1 


" 


= 











, / = 





1 





, A = 


1 





1 





















1 







1 


1 



and A 7 = I. 

We may replace the first four subspaces and the last subspace (which together contain 35 nonzero 
points) by ten two-dimensional subspaces with five single points left over. One choice for the ten 
two-dimensional subspaces is shown in Table [|. This gives a mixed spread of strength 2 consisting 
of four 3-dimensional subspaces, ten two-dimensional subspace and five points, and so by Theorem 



K)| corresponds to an 0,4(64, 2 5 4 1U 8 4 ). Finally, Lemma |l2| produces an 0,4(64, 2 5 4 17 8 lN 



Remark. Geometric orthogonal arrays with parameter sets of the form (64, . . . 8 k . . .) involve 
selecting k disjoint (except for the zero vector) copies of GF(2) 3 inside GF(2) G . It is simpler to 
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Table 4: Ten pairs of columns each spanning a two-dimensional subspace of GF(2) e . 



01 


01 





00 


01 


01 


00 


01 


10 


01 


00 


01 


01 


00 


01 


1 1 


01 


00 


01 


1 1 


01 


00 


10 


01 


00 


1 1 


00 


00 


01 


01 


00 


1 1 





00 


10 


01 


10 


10 


01 


10 


10 


01 


1 1 


00 


10 


01 


00 


1 1 


1 


00 


1 1 


10 


1 


10 


10 


01 


00 


00 


00 


1 



work projectively, and then we must choose k disjoint copies of PG(2, 2) inside PG(5, 2). Equation 
( | I ) then gives a decomposition of PG(5,2) into nine copies of PG(2,2). 

With the help of Magma (Bosma and Cannon, 1995; Bosma, Cannon and Mathews, 1994; 
Bosma, Cannon and Playoust, 1997), we showed that if 1 < k < 4 there is a unique way to choose 
k disjoint FG(2, 2)'s in PG(5,2), and these are equivalent to a subset of (|25|) . For k = 5, there 
are precisely two ways, one of which is equivalent to a subset of ( pq ) while the other contains no 
PG(2, 2) in its complement, and so cannot be extended to k = 6. For k = 6, . . . , 9, there is again a 
unique way to choose k disjoint planes. In particular the decomposition into nine planes shown in 
(|25| ) is also unique. 

An example of a maximal set of five PG(2, 2)'s in PG(5,2) is shown in Table [|. This corre- 
sponds to a geometric OA(64, 8 5 ) that cannot be extended to a geometric OA(QA, 8 6 ). It would be 
interesting to determine if it can be extended to any OA(6A, 8 6 ). 

Table 5: A set of five disjoint PG(2,2)'s in PG(5,2) that is not contained in a set of six. Each 
triple of columns spans one of the subspaces. 



100 


000 


100 


100 


100 


010 


000 


010 


010 


010 


001 


000 


001 


001 


001 


000 


100 


100 


010 


111 


000 


010 


010 


101 


110 


000 


001 


001 


011 


100 



7. The existence of orthogonal arrays with certain parameter sets 

In this section we prove Theorem [|, the lower bound in Theorem [2L and also give some other 
conditions which are sufficient to guarantee that a parameter set can be realized by an orthogonal 
array. 
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Lemma 14. Suppose N = p m is a power of a prime and (N,s kl S2 2 ■ ■ ■) is a parameter set with 
k = Yjiki factors. If k < pL( m+1 )/ 2 J + 1 then this parameter set is realized by a geometric orthogonal 
array. 

Proof. Suppose first that the parameter set contains a factor with s = p n > y/N levels. If m is 
even then a geometric OA(p m , (p m ~ n ) pn (p 71 ) 1 ) exists by Section 5, and p m ~ n is the largest number 
of levels other factors can have if there is an s-level factor. Since there are p n factors with p m ~ n 
levels, the existence of any array with one s-level factor and at most p m l 2 factors with < p m ~ n 
levels follows immediately. The case that m is odd follows similarly. 

We now assume that all Sj < y/N. If m is even, N = p 2r , then a geometric 

OA(p 2r ,(p r r r+1 ) (26) 

exists by Section 5. Any parameter set with all S{ < p r and k < p r + 1 is dominated by (Q]) and 
so is realized. If m is odd, N = p 2r+1 , then a geometric 

OA(p 2r+ \(p r ) pr+1+1 ) (27) 

also exists by Section 5. Any parameter set with all s» < p r and k < p r+1 + 1 is dominated by ([27]) 
and so is also realized. ■ 

Since the number of factors in a parameter set is less than or equal to the number of degrees of 
freedom (Eh, Lemma |l4| immediately implies that any parameter set (N = p" 1 ^ 1 ^ _ _ ^ w ^h at 
most pL(^+i)/2J + i de 

grees of freedom is realized by an orthogonal array. However, Theorem H is 

much stronger. 

Proof of Theorem We will show that any parameter set (p m , p kl (p 2 ) k2 (p^) ks . . .) satisfying 

J>(^ - 1) < P 3m/4 (28) 

is realized by a geometric orthogonal array, where p is any prime. To simplify the notation we 
assume m = 4r is a multiple of 4. The arguments in the other three cases require only minor 
modifications and are left to the reader. 
From ( p8[ ) we have 

k r+1 + k r+2 H + kur-i < p 2r + 1 , (29) 

and so by Lemma |l|| a geometric 

OA{p Ar , ^+1)^+1(^+2)^+2 _ _ _ ( p 4r-l)fc 4 r-i) 
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exists. 

We now proceed by induction. Let H n be the hypothesis that every parameter set 

(p 4r ,p bl (p 2 ) b2 . . . {p r ) br (p r+1 ) kr+1 . . . (p 4 *"- 1 )* 4 *- 1 ) (30) 
with hi < k{ for 1 < i < r and 

h + b 2 H \-b r = n 

can be realized by a geometric orthogonal array constructed using disjoint subspaces of PG{Ar — 
l,p). We have shown that Hq holds. Suppose H n holds with 

n < ki + k2 + • • • + k r . 

We will show that we can increase b r by 1 and still realize the parameter set, thus establishing 

To show this, note that the number of projectively distinct nonzero points in all the subspaces 
in ( |30"| ) is at most 

En — 1 p 
h 1 -— - 1 < - - 1 . 
p — 1 p — 1 

i=l 

However, by Theorem 1 of (Thas, 1995), Section 7, if a subset of PG(4r — l,p) contains fewer than 

(p 3r+1 - l)/(p- 1) points , 
there is a subspace PG{r — l,p) disjoint from it. Since 

1 < 



p — 1 p — 1 

such a subspace exists and we can use it to augment b r by 1. 

By induction, we can realize the parameter set (p 4r , p kl (p 2 ) k2 . . . (p 4r_1 ) fc 4r-i^ as required. ■ 

Proof of lower bound of Theorem First suppose N = 2 2r . From Theorem ^, every 
parameter set 

(2 2r 2 fcl 4 fc2 8 fc3 (2 r ) fcr ) 

with at most 2 3r / 2 degrees of freedom can be realized. The lower bound of Theorem |2| now follows 
in the same way that we proved the upper bound in Section 4. If N = 2 2r+1 we use the previous 
case together with T(N) > T(N/2). m 



It would be nice to have analogues of Theorem [1 and Lemma 14 when iV is not a prime power! 
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